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ABSTRACT 

The magnetic network which consists of vertical flux tubes located in intergranular 
lanes is dominated by Hall drift in the photosphere-lower chromosphere region (< 
1 Mm) . In the internetwork regions with weak magnetic field, Hall drift dominates 
above 0.25 Mm in the photosphere and below 2.5 Mm in the chromosphere. Although 
Hall drift does not cause any dissipation in the ambient plasma, it can destabilise 
the flux tubes and magnetic elements in the presence of azimuthal shear flow. The 
physical mechanism of this instability is quite simple: the shear flow twists the radial 
magnetic field and generates azimuthal field; torsional oscillations of the azimuthal 
field in turn generates the radial field completing feedback loop. The maximum growth 
rate of Hall instability is proportional to the absolute value of the shear gradient and is 
dependent on the ambient diffusivity. The diffusivity also determines the most unstable 
wavelength which is smaller for weaker fields. 

We apply the result of local stability analysis to the network and internetwork 
magnetic elements and show that the maximum growth rate for kilogauss field oc- 
curs around 0.5 Mm and decreases with increasing altitude. However, for 120 G field, 
the maximum growth rate remains almost constant in the entire photosphere - lower 
chromosphere except in a small region closer to the surface. For shear flow gradi- 
ent ^ 0.1 s~^, the Hall growth time is about 1 minute near photospheric footpoint. 
Therefore, network and internetwork regions with intense field in the presence of shear 
flow are likely to be unstable in the photosphere. The weak field internetwork regions 
could be unstable in the entire photosphere-chromosphere. Thus the Hall instability 
can play an important role in generating low frequency turbulence which can heat the 
chromosphere. 
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1 INTRODUCTION 

Although our understanding of solar magnetism remains 
incomplete with each observational improvement revealing 
new magnetic structures, the following canonical picture re- 
mains unaltered. The photosphere is threaded by strong ver- 
tical magnetic field concentrated in individual elements or 
flux tubes at intergranular boundaries (Simon & Leighton 
1964). The field strength within individual magnetic ele- 
ments are > 1 — 2 kG and their typical size is ~ 100 — 200 km 
at the footpoints. Clear spatial correlations between Call K 
line, UV intensity and heating in the chromosphere maps out 
the underlying photospheric fields quite well (Hasan 2009; 
Loukitcheva et al. 2009). The magnetic network shows up as 
a collection of Call bright points in the H and K lines and its 
appearance is taken as a proxy for foot point motion of the 
flux tubes (Hasan et al. 2000; Hasan 2008). The lifetime of 
bright points varies between ~ 200 to 400 s (Krishnakumar 



& Venkatakrishnan 1999). Isolated intense field concentra- 
tions that produce bright points have also been detected in 
the internetwork region (de Wijn et al. 2009; Almeida et 
al. 2010). Outside the concentration, non- vertical (possibly 
isotropic, Almeida & Gonzalez, 2011) weaker field also ex- 
ists in the interior of supergranules (Lites et al. 2008). The 
flux tubes fan out from the photosphere (with filling factor 
< 1%) to chromosphere where individual tubes merge and 
form a canopy (with filling factor ~ 100%). 

The search for a viable mechanism that heats the chro- 
mosphere has been underway for the past several decades. 
The original idea of acoustic heating due to photospheric 
granular motion (Biermann 1946; Schwarzschild 1948; Stein 
1967; Carlsson & Stein 1992; Lites et al. 1993; Kalkofen 
2007) implies that such a mechanism could be important 
in regions where the magnetic field is dynamically unimpor- 
tant. Observations by NASA's Transition Region And Coro- 
nal Explorer (TRACE) satellite revealed that there is hardly 
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any acoustic energy flux (about 90 % deficit) to heat the non- 
magnetic, quiet chromosphere (Fossum & Carlsson 2005a,b, 
2006). Subsequent numerical simulations suggested that this 
huge discrepancy in the energy flux is possibly due to the 
limited spatial resolution of TRACE (Wedemeyer-Bohm et 
al. 2004; Cuntz et al. 2007). Recent high-resolution observa- 
tions suggest that the total acoustic power is at least twice 
as large as the largest power observed before (Gonzalez et al. 
2010). However, the jury is out on various acoustic heating 
models of the chromosphere. 

As intense heating in the chromosphere 
(~ 10^ ergs cm~^ s~^) is directly correlated with the 
strongest magnetic fleld concentrations in the network and 
internetwork regions (Schrijver et al. 1989; Loukitcheva et 
al. 2009), the magnetic fleld must play important role in 
heating the chromospheric plasma. For example, resistive 
heating of the plasma due to Pedersen diffusion is one such 
possibility (Goodman 2004). However, recently it has been 
pointed out that resistive heating may not be dominant 
in the chromosphere (Socas-Navarro 2005). Drawing upon 
the similarity between Earth's ionosphere and the Solar 
atmosphere, the Farley-Buneman instability, which is 
possibly responsible for the equatorial electrojet and can 
be triggered by cross field motion of the partially ionized 
medium, has been proposed as another plausible heating 
mechanism of the chromosphere (Fontenla et al. 2008). It is 
quite likely that the fiow threshold, > 2km/s, required to 
drive this electrostatic instability may not be achievable in 
the chromosphere (Gogoberidze et al. 2009). Therefore, the 
problem of chromospheric heating remains unsolved. 

It is well known that Hall diffusion can not directly 
heat the plasma and thus has been overlooked as a possible 
cause of network-internetwork heating. However, Hall diffu- 
sion in the presence of velocity shear can destabilise MHD 
waves (Riidiger & Kitchatinov 2005; Kunz 2008; Bejarano 
et al. 2010; Wardle & Salmeron 2012) and thus can indi- 
rectly heat the medium by driving turbulence. The crucial 
elements required to excite Hall instability are (a) the pres- 
ence of inhomogeneous azimuthal flow, and, (b) Hall drift of 
the magnetic fleld. It is quite likely that both these condi- 
tions are met in the photosphere-chromosphere. 

Rotation has been invoked in the past to explain the 
stability of flux tubes (Schiissler 1984). Models of spicules 
also use the concept of rotating flux tubes (Kudoh & Shi- 
bata 1997). The swirling macrospicules, manifested through 

o 

blue and red — shifted emission in the OV line (629.73 A) 
at the transition region indicate rotational velocities ~ 
20 - 30km/s (Pike & Mason 1998). The recent discovery of 
convectively driven localised vortex-type motions (diameters 
< 500 km; lifetimes 3 — 7min.) which occur in intergranu- 
lar magnetic concentrations in the photosphere (Bonet et al. 
2008) has raised the possibility that it may play important 
role in the energy transfer from photosphere to corona. The 
analysis of SUNRISE observations (Bonet et al. 2010) sug- 
gests the typical lifetime of the vortices is between 5 and 20 
minutes with a mean ~ 7.9 ±3.2 minutes. Wedemeyer-Bohm 
& Voort (2009) analysed the ionized calcium line and found 
the presence of rotation in the chromosphere which they in- 
ferred to be a manifestation of the twisting and braiding of 
magnetic footpoint in the photosphere. 

Numerical simulation of near-surface solar convection 




Height [Mm] 

Figure 1. The ratio of Pedersen (rjp) and Hall (i^h) diffusion is 
shown in the above figure. The diffusivities have been calculated 
for magnetic field profile B = Bo(n„/no)'' '' where Bq = 1.2kG , 
(bold line), and 120 G , (dotted line), n„ is neutral number density 
and no is neutral number density at the solar surface, h = 0, taken 
from model C, Vernazza et al. (1981). 

displays turbulent vortex flows at intergranular lanes (Zirker 
1993; Stein & Nordlund 1998). The observation of granu- 
lation in the very quiet solar region in the disc centre sug- 
gests the presence of vortex tubes with typical mean radius 
~ 150 km and rotation period ~ 30 s considerably shorter 
than the lifetime of a typical vortex tube (~ few minutes) 
(Steiner et al. 2010). Vorticity generation near boundaries of 
granules has also been seen in numerical simulations of the 
photosphere (Stein & Nordlund 1998; Muthsam et al. 2010). 
Small-scale, intergranular vortices are formed due to inter- 
action of the photospheric plasma with ambient magnetic 
fleld (Moll et al. 2011; Shelyag et al. 2011). Recent high res- 
olution simulations of a realistic photosphere shows that the 
Hall effect generates out-of-plane velocity flelds with max- 
imum speed ~ O.lkm/s at interface layers between weakly 
magnetized light bridges and neighbouring strong field um- 
bral regions (Cheung & Cameron 2012). To summarise, both 
observation and numerical simulation point to the frequent 
occurrence of shear flows in the solar photosphere. There- 
fore, both of the above mentioned conditions necessary for 
the onset of Hall instability are likely to be met in the 
photosphere-chromosphere. 

In Fig. (1) we show the ratio of Pedersen, and. Hall 
diffusivities, rjp/rjH- Although both ion and neutral num- 
ber densities decrease with increasing height, the ratio B /rii 
where rii is ion number density do not always decrease and 
this is reflected in jagged diffusivity proflles.^ Fig. 1 sug- 
gest that Hall dominates Pedersen diffusion below ~ 1 Mm 
in the network (~ kG) region. In the internetwork region, 
where the fleld is ~ 100 G, Hall diffusion dominates in the 
entire photosphere above 0.25 Mm and the chromosphere. 
Thus, conditions (a) and (b) required for the onset of Hall 

^ The detailed calculation of various diffusivities for the 
photosphere-chromosphere region will be given in Pandey &; War- 
dle (2012). 
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instability are easily met in the solar atmosphere. Clearly, 
Hall diffusion in the presence of shear flow can destabilise 
the flux tubes and magnetic elements and drive turbulence 
in the medium which could facilitate the heating. 

In the present paper, we investigate the local stability of 
the network and internetwork magnetic elements by assum- 
ing that magnetic fleld is immersed in the highly diflxisive 
plasma medium. Since the solar atmosphere is highly strat- 
ified, validity of present local analysis is restricted to short 
(with respect to the scale height) wavelength fluctuations. 
Further, since flux tube will be approximated by a planar 
geometry, our analysis is valid only for wavelengths much 
smaller than the tube radius. 



+ V ■ (pv) = 0. 



(1) 



dt c 
where J = eue {vi — Ve) is the current density, B is the 
magnetic field and P — + Pi + Pn is the total pressure. 

The induction equation is 



dB 
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where b = B/B is the unit vector along magnetic field. 

The expression for Hall, rjH, ambipolar, rjA and. Ohm, rj 
diffusivitics arc (Pandey & Wardlc 2008) 
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where D = pn/ p , va ~ B / ^/TiTp is Alfvon velocity, /3e = 
Wee/ Ven is the elcctrou Hall parameter, a ratio of the electron 
cyclotron to the electron-neutral collision frequencies and 
the Hall frequency u)h is 



Pi 

Wif = — We 

P 



J^e We 



(5) 



with Xe = Ue/rin. 

The concept of frozen-in flux is quite useful in ideal 
MHD and thus it would be profitable to cast the induction 
Eq. (3) in similar form. Such a formulation allows us to 
visualise the magnetic field as a real physical entity that is 
drifting through the fluid with a given velocity. Lets consider 
magnetic flux ^ = J B ■ dS through surface S encircled by 
an arbitrary closed contour C in the plasma moving with 
velocity u. The time rate of change of the magnetic flux is 
given as 

'^^ f ^-^8+ (f B {uxdl) 



dt 



E+-UX B 



dl 



(6) 



2 BASIC MODEL 

The partially ionized photosphere-chromosphere plasma 
consist almost entirely of electrons, protons, singly ionized 
metallic ions, H, He I, He H, and He HL Wc shall assume 
that the plasma consists only of electrons, hydrogen ions and 
neutrals and neglect the distinction between the metallic and 
hydrogen ions. The single-fluid dcscriptiou of partially ion- 
ized plasma is well known for past several decades (Cowling 
1957; Braginskii 1965; Mitchner& Kruger 1973). However, 
only recently Pandey & Wardle (2008) critically analysed 
the non-ideal MHD effects in a partially ionised plasma by 
developing a unified singlc-fiuid framework for the plasmas 
of arbitrary ionization. Therefore, we shall utilise single fiuid 
formulation given by Pandey & Wardle (2008). Note that 
plasma quasineutrality is implied in the single fluid formu- 
lation. 

The continuity equation is 

dp 
dt 

Here p !^ pi + pn is bulk fluid density and v = {pi Vi + 

pn Vn)/p is bulk velocity with pi, Vi and p„, Vn representing 
the mass density and bulk velocities of the ion and neutral 
fluids respectively. The momentum equation is 

/^ = -VP+^, (2) ^ 



where dS is the surface element and dl is the line element 
along C. The flux $ through the closed contour C is con- 
served only when the integrand in the last expression of 
Eq. (6) vanishes, i.e. cE + uxB = 0. 

We note that the induction Eq. (3) is derived from gen- 
eralized Ohm's law 



—E' = 'nJ«+riHJxb + ripJ±. 

An 



(7) 



where electric field is written in the neutral frame and paral- 
lel and perpendicular components of the current J refers to 
the orientation with respect to the ambient magnetic field 



( J • b) b , Jj_ = J - J|| 



(8) 



The Pedersen diffusion r)p is related to the ambipolar 
diffusion as 

rjA = VP - V ■ (9) 

Thus defining magnetic drift velocity as (Wardle & Salmeron 
2012) 

(V xB)^xb (V X B)^ 
VB = rjp riH ^ , (10) 

the induction Eq. (3) can be explicitly written in terms of 
fluid and field velocities as 

dB 



dt 



Vx 



4 TT 77 

{v + vb)xB '- J|| 



(11) 



We may infer from cE± + {v + vb)xB = that when 
each point of the boundary C moves with the velocity 



■V + Vb, 



(12) 



the time rate of change of flux in partially ionized medium 
becomes 

'dt 



(VxS)ii • dl. 



(13) 



Wc see from Eq. (13) that the presence of parallel current 
is responsible for the flux non-conservation although rate of 
flux decay is directly proportional to Ohm diffusion. In the 
absence of parallel current, non- ideal MHD effects only re- 
distribute magnetic flux in the medium leaving total flux un- 
changed. At first sight this result appears paradoxical since 
non-ideal MHD effects, particularly, Pedersen (Ohm -|- Am- 
bipolar) diffusion, owing to its dissipative nature, must af- 
fect magnetic flux. Indeed energy dissipation is inevitable 
in such collision dominated weakly ionized plasmas. How- 
ever, energy loss occurs due to redistribution / relaxation of 
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Figure 2. The magnetic Reynolds numbers Rm, Am and Hm 
pertaining to Ohm, ambipolar and Hall diflfusion is shown against 
height for photosphere-chromosphere plasma parameters of model 
C, VAL81. 



the magnetic field and not due to flux annihilation (Parker 
1963). Therefore, if parallel current is absent in the medium, 
the total flux is conserved. 

Although for the past two decades, the dynamics of the 
weakly ionized-weakly magnetized photosphere and weakly 
ionized-strongly magnetized chromosphere has been inves- 
tigated in the framework of non-ideal MHD [Arber et al. 
(2007), Goodman (2011) and references therein] it is per- 
tinent to estimate relative importance of various terms in 
the induction equation (3) . They are on the order of ct B : 
vB/L : -nB/L^ : t)hB/L^ : r]AB/L^ . Here we have re- 
placed dt by a which is dynamical frequency set by the sig- 
nal speed and gradient length scale L (Pandey & Wardle 
2008) . The ratio of convective and diffusive terms in the in- 
duction Eq. (3) is known as Magnetic Reynolds number Rm 
in MHD (Goedbloed & Poedts 2004). We shall generalise 
this definition of magnetic Reynolds number and, introduce 
three different magnetic Reynolds numbers corresponding to 
Ohm, ambipolar and Hall diffusion as 

Rm = , Am = , Hm = . 

The ambipolar and Hall Reynolds numbers Am and Hm 
can be expressed in terms of magnetic Reynolds number 
Rm with the help of ion and electron Hall parameters. 
Thus defining ion Hall parameter Pi — tod/vin, where 
ujci = eB/rriiC is the ion cyclotron frequency, ambipolar 
and Hall Reynolds numbers can be written as 



(14) 



Am = Rm/ {|3^ P^) , Hm = Rm//3e 



(15) 



Clearly, both ambipolar and Hall Reynolds numbers are 
functions of plasma magnetisation, that is plasma Hall pa- 
rameter /?e and Pi (Pandey et al. 2008) and therefore, their 
values in weakly ionized and weakly magnetized photosphere 
and weakly ionized and strongly magnetized chromosphere 
widely differ. 

In Fig. (2) we plot Rm, Hm and Am against height for 
a kG [Fig. 2(a)] and 0.1 kG [Fig. 2(b)] fields. As has been 
mentioned earlier, we have utilised model C, Vernazza et al. 
(1981) (hereafter VAL81) to calculate various diffusivities. 



Here we assume L — 100km and v = va- We see that for 
both kG and sub kG fields, magnetic Reynolds number Rm 
is very large except close to the surface. Therefore, barring 
the solar surface. Ohm diffusion can be neglected in the en- 
tire photosphere-chromosphere. However, unlike Rm, Hall 
and ambipolar Reynolds number are not very large for a kG 
field in most of the photosphere-chromosphere. In fact in 
the large part of photosphere-chromosphere (> 0.25 Mm), 
Hm takes values between lO'^ and 10^ . Similarly, Am varies 
between 1 and 10'^ in this interval. Clearly, these quanti- 
ties are several oreder of magnitude smaller than Rm. The 
present estimate of Hm and Am is broadly in agreement 
with Khomenko & CoUados (2012); Sykora et al. (2012). 
For sub-Alfvenic motion (Hiraki et al. 2010), both Hm and 
Am can become of the order of unity or less. For 0.1 kG 
field [Fig. 2(b)], advection term in induction Eq. (3) is more 
dominant than Hall and ambipolar diffusion in comparison 
with the previous case. However, at smaller scales or, sub- 
Alfven speed, even for a weak field, Hall and ambipolar dif- 
fusion becomes important. 

Recent radiative MHD simulation of weakly ionised so- 
lar atmosphere shows that Hall term is largest in the lower 
and middle chromosphere and in the corona. In the inter- 
granular lanes in the photosphere, the Hall term is the most 
important diffusion term whereas ambipolar diffusion is im- 
portant in the region from the upper-photosphere to the up- 
per chromosphere (Sykora et al. 2012). In the chromosphere, 
ambipolar diffusion dominates almost everywhere except in 
the lower chromosphere in shock fronts. To conclude, both 
Hall and ambipolar diffusion can compete with the advection 
term in the induction equation and thus, must be retained 
while investigating the weakly ionized solar atmosphere. 



3 STABILITY ANALYSIS 

As magnetic elements in the solar atmosphere are concen- 
trated into discrete structures, they are often modelled as 
flux tubes or, flux slabs, implying that the field either posses 
cylindrical or translational symmetry. Further, both in ac- 
tive and quiet phases, magnetic elements are highly dynamic 
consisting of numerous flows with different spatial and tem- 
poral time scales. Recent two dimensional simulations of um- 
bral magneto convection suggest that the dynamical scale 
over which Hall effect can generate magnetic and velocity 
fields is much faster [~ 10-20 km spatial scale and ~ 300 s 
temporal scale, Cheung & Cameron (2012)] than the spa- 
tial and temporal scales of the flux tubes (> few hundred 
km and ~ few days). In fact, the numerical results are eas- 
ily scalable to 2 km with temporal scale ~ 2 s. It should be 
pointed out that present observations can not resolve the 
shear flow scale since best achievable resolution is of Hinode 
is > 100 - 200 km. 

As has been noted above, the spatial scale over which 
flow and field generation occurs is much smaller than the 
typical tube diameter. Thus we shall approximate the cylin- 
drical tube by a planar sheet and work in the Cartesian co- 
ordinates where x ,y ,z represent the local radial, azimuthal 
and vertical directions. We shall assume an initial homo- 
geneous state with azimuthal shear flow v = vo'xy and 
uniform vertical field, B — B z. We note that coUisional 
dissipation will cause the loss of energy in the lower solar 
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atmosphere and in general a proper energy equation should 
be used for a more realistic modelling of the physical pro- 
cesses (Khomenko &: Collados 2012). However, in order to 
keep the description simple, and to understand the under- 
lying physical processes we shall assume that the fluid is 
incompressible. 

We shall assume that the solutions of the linearised 
equations is of the form cxp (a t + ikz). This simple form 
allows us to split the linear system of equations in two sub- 
systems: one subsystem corresponding to sound waves prop- 
agating along the magnetic field and the interesting subsys- 
tem for which z — components of the perturbed velocity, 
magnetic fleld, current density and electric field all vanish. 

The X and y components of the momentum equation 
yield 



Sv 

VA 



i k 



a 
a a 



5B 
B 



The magnetic drift velocity, Eq. (10) becomes 
Svb 

VA 



■ i k 



VP Vh 

-m VP 



SB 

B 



(16) 



(17) 



Since 5J\\ = 0, x and y components of the induction Eq. (11) 

can be written as 



a 
a a 



6B .rfSv Svb 
— =ik — + 

B \VA VA 



(18) 



In the above Eqns. (16)-(18), wavenumber k, frequency a 

and diffusivity rj are normalised: k = kvA/\v'o\, a = (y/\vo\ 
and 77P = r]p \v'o\/va, 7)h = rjH \v'o\/va ■ Here va = B/y/AiTp 
is Alfvenspcod and a — —v'q/ \vq\ = ±1. Making use of 
Eqs. (16)-(17), Eq. (18) can be written as 



+ k' 





(T 
f/p 

-m 



VP 



SB 

B 



= 0. 



(19) 



Since the dynamical time of interest here is set by shear 
flow gradient i.e. cr ~ \v'o\ or, a ■-^ 1 three terms in the 
preceding equation is of the order of 1 , k^ , k^ fj±, where 
fj± = v^H+^p- Thus, Eq. (19) can be analysed in three 
different limits (Wardle & Salmeron 2012). 

I. Ideal MHD: ^ Magnetic diffusion is negligible in com- 
parison with the fluid advection, i.e. the field is frozen in 
the fluid; [Sv/va ~ SB/B]. In this limit fc^ ~ 1 > fe^j^x, 
i.e. fj± <^ 1 and the last term in Eq. (19) can be neglected. 
Since fc'^ ~ 1 implies kvA ~ |fol for va = 5km/s and vor- 
tex flow |t;o| ~ 0.1/s (Stein & Nordlund 1998), this gives 
A ~ 300 km. The pressure scale height H in the photosphere 
is ~ 150 km (Spruit 1976). Thus A > and the ideal MHD 
limit is not applicable to the thin flux tubes. 

II. Cyclotron limit: Magnetic diffusion is comparable to 
the fluid advection; [5vjvA ^ SB/B]. In this case field dif- 
fusion balances fluid convection, i.e. fc^ v± ~ fc^ > 1. This 
is the low frequency limit since ct ~ 1 and first term in 
Eq. (19) can be neglected. The low frequency, short wave- 
length dressed ion-cyclotron wave with frequency ujc = ijJh 
is the normal mode of the system (Pandey & Wardle 2008). 

^ Here ideal MHD of weakly ionized plasma pertains to the neu- 
tral fluid which carries the inertia of the medium. 



The requirement k^ ii± ^ k^ ^ 1 is never fulfilled in the 
solar atmosphere since v± = v± I'^ol/'^i ~ 10~^ — 10~* -C 1. 
Therefore, this limit is not valid in the solar atmosphere. 

HI. Highly diffusive limit: Magnetic diffusion over- 
whelms fluid advection; [Sv/va <S SB/B]. In this case 
V± ~ 1 ^ k^, i.e. fj± ^ 1 and k^ <C 1. The evolution 
of magnetic field which is kinematic in nature in this case is 
solely determined by diffusion. Only the first and last terms 
in Eq. (19) are retained. For typical values of v± (Table 1), 
this limit gives A < 6 km which fits within the pressure scale 
height. Therefore, we shall work in the high frequency limit 
and neglect the advection term in Eq. (19). 

In the highly diffusive limit, we get the following dis- 
persion relation 



a+Cia + Co = 0, 
where 

Ci = 2 k'' Tjp , Co = fe^ [fc^ ri'i -I- Tin] 



(20) 



(21) 



This dispersion relation can as well be written in non- 
dimensional form a'^ -\- Ci a -\- Co = , with coefficients 



Ci = 2 fe^ 7?P , Co = k^ 



k r/x - a rju 



(22) 



In the absence of shear flow and VPj fhe right circularly 
polarised whistler 



ojw = {kvAf /ojh = k'^ vh ■ 



(23) 



is the normal mode of the system (Pandey & Wardle 2008). 

The necessary condition for Hall instability, Co < 0, 
becomes 



-v'o Vh > k'^ v±- 



(24) 



Note that the sign of Hall diffusion (vh) depends on the 

orientation of the vertical magnetic field with respect to 
the shear flow gradient (Wardle 1999). As a result onset 
of the Hall instability depends not only on the sign of the 
velocity gradient but also on the orientation of the vertical 
magnetic field. Thus when the shear gradient is such that 
—v'o > 0, a positive rjH (which can be guaranteed if the 
magnetic field is parallel to the shear gradient, —v'q) allows 
the above condition to be easily satisfied. When the mag- 
netic field is antiparallel to the shear gradient v'q > 0, above 
condition can as well be satisfied. Therefore, the orienta- 
tion of the magnetic field together with the shear gradient 
plays crucial role in the Hall instability. This instability is 
the planar version of well known Hall-modified magnetoro- 
taional instability in accretion discs where the differential 
rotation of the disc [v{x) = xQ{x) where n is the orbital 
frequency] causes shear in the disc (Wardle 1999; Wardle & 
Salmeron 2012). For definiteness in the subsequent analysis 
we shall assume a = 1. 

Eq. (24) suggests that the fiuctuations of wavelength 



A > 



2tv 



vl 

VH + — 
VH 



1/2 



(25) 



will become unstable. Flow gradients with either sign will 
give real A provided magnetic field has the right orientation. 
The wavelength > 60 km becomes Hall unstable for r]H ~ 
10^° cmVs and vp ~ Vh (Table 1). 

The growth rate of Hall instability is 
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Figure 3. The growth rate of Hall instability a = <y/\v'g \ against 



k imy/^ = k (riHl 
VP/VH- 



,1/2 



is shown for different values of 



(26) 



In Fig. (3) Eq. (26) is plotted for various values of 
Vp/Vh and a — 1. We see that with increasing rjp the 
growth rate decreases before completely disappearing when 
rjp S> rjH, that is when the available free energy is com- 
pletely dissipated by Pedersen diffusion. 

The physical picture of Hall instability is quite simple. 
The shear flow generates SBy from SBx and Hall diffusion 
generates SB^ from SBy. This is how Hall diffusion in tan- 
dem with shear flow destabilises the flux tube. This also ex- 
plains the dependence of growth rate on both shear gradient 
and whistler frequencies. 

We can find the maximum growth rate of Hall in- 
stability by recasting the dispersion relation, Eq. (22) in 



4 a c = 0. Here 



-2 7 



2fipa — fjn, c 



The maximum growth rate of the instability becomes 



VP + r]i_ 



2 



np/ilH + \/l + ril/rin 



(27) 



.(28) 



It is clear from the preceding equation that the maximum 
growth rate is proportional to the shear gradient and to lead- 
ing order, is inversely proportional to 1 -I- rfp/riH- Thus, in 
conformity with Fig. (3), the maximum growth rate of the 
instability corresponds to rjp/riH = and with increasing 
Vp/Vh, the instability grows at smaller rates. The wavenum- 
ber corresponding to maximum growth rate is 



'2il± (Vp + V±) 



(29) 



Thus the maximum growth rate is ctq = fco??x. Note that 
fco oc 1/ ^y^r|H ~ 1 and hence for a given rip /'qu , a weaker 
field implies that smaller wavelengths are Hall unstable. 



4 APPLICATION TO THE SOLAR 
ATMOSPHERE 

Intense magnetic fields are organised as mainly vertical fiux 
tubes and are thought to attain steady state possibly due 
to the pressure balance with the ambient medium. Cylin- 
drical coordinates are best suited to describe these struc- 
tures. However, in order to develop proper understanding of 
physics behind the Hall diffusion driven shear flow instabil- 
ity, we have approximated cylindrical ffux tubes locally by 
planar sheets in the present work where x and y coordinate 
locally correspond to the radial and azimuthal directions re- 
spectively on a cylinder. 

For a kiloGauss field. Hall diffusion dominates photo- 
sphere and lower chromosphere between and 1 Mm (see 
Fig. 1). Thus, in active and quiet solar regions closer to 
the footpoint {h = 0) where kG field may be present. Hall 
instability will have the maximum growth rate which is 
on the order of the shear gradient [r^p = in Fig. (3)]. 
With increasing height, for example at ~ 1.5 Mm, when 
Tjp/riH ~ 4 (Fig. 1), the Hall instability growth rate is one 
third of shear gradient. At ft ~ 2 Mm, when rjp/rjH = 1-75 
the growth rate reduces to ~ 0.35 |wo|. Therefore, with in- 
creasing height (> IMm), when Pedersen diffusion domi- 
nates Hall diffusion. Hall instability grows at considerably 
smaller rate than below 1 Mm. 

In the internetwork regions, when Bo ~ 100 Gauss, 
Hall diffusion dominates Pedersen in the entire photosphere- 
chromosphere except for a very small region closer to the 
solar surface (< 0.25 Mm) (dotted curve in Fig. 1). As a 
result, Hall instability will grow at a rate close to shear fre- 
quency (curves labelled and 0.2 in Fig. [3]), in the entire 
photosphere-chromosphere above 0.25 Mm. 

The maximum growth rate of Hall instability is inde- 
pendent of the magnetic field strength only when rjp /rjH <C 
1. This implies that for the kG network or internetwork field 
below 1 Mm or, weak internetwork fields above 0.25 Mm 
(where Hall is the dominant diffusion), Hall instability 
will not be affected by decreasing magnitude of the field. 
However, in the weak-field internetwork regions below < 
0.25 Mm, where rjp/rjH ^ 1 the maximum growth rate can 
be quenched by the presence of a strong field, since rjp ~ B^ 
and rjH ~ B, and thus the maximum growth rate is inversely 
proportional to the magnetic field. Therefore, in strong field 
regions, long wavelength fiuctuations are suppressed by Ped- 
ersen diffusion. 

In the present problem we have assumed that the par- 
tially ionized plasma is threaded by a uniform vertical mag- 
netic field. This is the simplest representation of magnetic 
fiux tubes in the solar atmosphere. We could directly ap- 
ply these results to fiux tubes provided the fiow velocity 
and shear scale were known. Therefore, we turn to recent 
observation and numerical simulations. The vortex fiows 
on the solar surface have been observed by various groups 
[Please see excellent review by Steiner & Rezaei (2012)]. 
Small whirlpools similar to terrestrial hurricanes (< 500 km) 
with mean lifetime (~ 5min.) were discovered by Bonet et 
al. (2008). Larger vortices with longer lifetime and possi- 
ble entanglement of magnetic field manifested as enhanced 
Call emission at vortex centre have also been observed 
(Attie et al. 2009). The typical vorticity of a vortex is 
which corresponds to rotation period ~ 35 
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Height [Mm] 

Figure 4. The maximum growth rate and most unstable 
wavenumber normalised to pressure scale height for 1.2 kG (bold 
line) and 120 G field (dotted line) against altitude measured from 
solar surface [h = 0) is shown in Fig. 4(a) and Fig. 4(b) respec- 
tively. 

minutes (Bonet et al. 2010). Thus it would appear that 
Hall instability does not have time to develop since the 
growth rate {\v'o\/2 = 3 x 10"^ s"!) is very small. However, 
above vorticity value is limited by the upper limit in the 
vorticity resolution (~ 4 x 10"^ s~^, Bonet et al. (2010)). 
The numerical simulation gives much higher vorticity value 
(~ 0.1 — 0.2 s~^) in the photosphere-lower chromosphere 
(Fig. 31, Stein & Nordlund (1998)). The growth rate corre- 
sponding to \v'o\ = 0.2 is one minute. 

For a kilo Gauss field at the footpoint of the flux tube 
{h = 0), using corresponding number densities of neutrals 
and ions (model C, VAL81), the Alfven speed becomes va ~ 
5 X lO'^ cm/s. In Table 1, we give values of Hall diffusion the 
photosphere-chromosphere region which has been used to 
estimate fluctuation wavelengths. Since rjH ~ B, the value 
of Hall diffusivity in the internetwork region with weak field 
~ 120 G will be an order of magnitude smaller than values 
given in the table. 

In Fig. [4(a)-(b)], we plot the maximum growth rate, 
Eq. (28) and corresponding wavelength Eq. (29) for |wo| — 
0.1 s~^ against altitude. As expected the maximum growth 
rate (~ \v'o\) for kilo Gauss field occurs around 0.5 Mm, 
where rjp/rjH ^ 1 (Fig. 1). With increasing altitude, since 
Pedersen becomes comparable to Hall diffusion, the growth 
rate tapers off and beyond ~ 1 Mm Hall instability may 
not be dynamically important. However, for a 120 G field, 
excluding a small region in the lower photosphere, the 
maximum growth rate remains robust, < \v'(^\ in entire 
photosphere-chromosphere. 

In Fig. 4(b) we plot the most unstable wavelength nor- 
malised against pressure scale height H which has been 
calculated using model C, VAL81. For both intense and 
weak fields, Ao < H and the fluctuation wavelength fits 
within the scale height. Therefore, the most unstable wave- 
length fits within the scale height in the entire photosphere- 
chromosphere. 

It would appear from Eq. (22) that the wavenumber 
spectrum has finite cutoff. However, cut-off wavelength can 
not be inferred from Eq. (22) as cr ~ 1 and k ~ in 
highly diffusive limit. In the cr — >■ limit, the neglected ad- 
vection term in the induction equation also becomes im- 



portant. Therefore, we can not infer cut-off wavelength in 
the highly diffusive limit. The spectrum of Hall instabil- 
ity depends on the ambient magnetic field strength, i.e. 
fco oc Xj^fB implying that the wave spectrum will have 
sharper peak in kG field regions. 

With increasing altitude vertical tubes expand and 
bend developing small azimuthal field components in the 
process. As will be shown elsewhere (Pandey & Wardle 
2012), growth of the Hall instability is the merely rescaled 
due to bending of the field. 



5 DISCUSSION AND SUMMARY 

Granular motions are responsible for the generation of low 
frequency waves in the predominantly neutral photosphere 
and lower chromosphere. In the presence of a shear flow gra- 
dient, tubes in both network and internetwork regions can 
become unstable due to Hall diffusion. Recall that the insta- 
bility growth rate is «o/2 which for iig = 0.1 s"'^ correspond 
to the e-fold time ~ 20 s. Since average period of vortical 
motion in the photosphere is ~ 300 s. Hall instability have 
sufficient time to develop in the vortices. The observed value 
of vorticity is an order of magnitude smaller (Bonet et al. 
2010) than numerical value suggesting that Hall instability 
may not have enough time to develop. However, observed 
small value of the vorticity is due to limited spatial-temporal 
resolution (Bonet et al. 2010). The simulation results indi- 
cate that highest vorticities occur at the smallest observable 
scale (Stein & Nordlund 1998). Therefore, we conclude that 
Hall instability in tandem with shear flow may destabilise 
the flux tubes. 

Can Hall instability excite turbulence and heat the 
medium? This question falls beyond the scope of present lin- 
ear analysis. Only numerical simulation can provide an an- 
swer to this question. Numerical simulations are approach- 
ing length scales associated with the required shear scale 
(MoU et al. 2011; Cheung & Cameron 2012; Sykora et al. 
2012). Therefore, we hope that simulation by various groups 
will be able to delineate the role of Hall instability in chro- 
mospheric heating. The present analysis provides plausible 
pathway to such a process. 
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